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We compute the finite temperature Casimir energy for massive scalar field with general curvature 
coupling subject to Dirichlet or Neumann boundary conditions on the walls of a closed cylinder with 
arbitrary cross section, located in a background spacetime of the form M'*^"'"'^ x A/"", where Af^^^^ 
is the (di + 1) -dimensional Minkowski spacetime and A/"" is an n-dimensional internal manifold. 
The Casimir energy is regularized using the criteria that it should vanish in the infinite mass limit. 
The Casimir force acting on a piston moving freely inside the closed cylinder is derived and it is 
shown that it is independent of the regularization procedure. By letting one of the chambers of 
the cylinder divided by the piston to be infinitely long, we obtain the Casimir force acting on two 
parallel plates embedded in the cylinder. It is shown that if both the plates assume Dirichlet or 
Neumann boundary conditions, the strength of the Casimir force is reduced by the increase in mass. 
Under certain conditions, the passage from massless to massive will change the nature of the force 
from long range to short range. Other properties of the Casimir force such as its sign, its behavior 
at low and high temperature, and its behavior at small and large plate separations, are found to 
be similar to the massless case. Explicit exact formulas and asymptotic behaviors of the Casimir 
force at different limits are derived. The Casimir force when one plate assumes Dirichlet boundary 
condition and one plate assumes Neumann boundary condition is also derived and shown to be 
repulsive. 
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I. INTRODUCTION 



In the endeavor to solve some fundamental problems in physics, such as the unification of fundamental forces and 
the dark energy and cosmological constant problem, it has been proposed that we should consider higher dimensional 
spacetime. For example, string theory [l[ predicts that we live in spacetime of ten or eleven dimensions, where the 
extra six or seven space dimensions are curled up to a tiny invisible compact manifold. Therefore, there is a strong 
motivation to study physics in spacetimes with extra dimensions. Since Casimir effect is a fundamental quantum effect, 
the influence of extra dimensions on Casimir effect becomes an important issue. The Casimir effect in spacetime with 
extra dimensions was considered in the works 3, 1^, B [ a jjO, 11, l3. [Tsl . [T3. [Tsl. [l^ [iTj . The role of Casimir 
effect in stabilizing extra dimensions was discussed in [1^, [l^, [20, [2l|, (22, [23, [2J] ■ The possible role of Casimir effect 
as cosmological constant responsible for the observed dark energy was considered in O, [H, [H [13, [H [H, [11 . In 
[3) [E]j the Casimir effect for massless scalar field in Kaluza-Klein spacetime of the form M^+^ x A/"", where M^+^ 
is the (3 -|- l)-dimensional Minkowski spacetime and A/"" is an n-dimensional internal manifold, was considered in 
the piston setting. In [iJi [IB| , we extended the results of [1, [^ and considered the finite temperature correction to 
the Casimir effect. In the work [l6|, the Casimir effect due to massive scalar field with general curvature coupling 
constant subject to Robin boundary conditions in spacetime of the form M'^^'^^ x A^", where M'''^^^ is the [di + 1)- 
dimensional Minskowski spacetime, was considered. In this paper, we generalize the work of (l6j by taking into 
account the temperature correction. However, we restrict ourselves to only consider either Dirichlet or Neumann 
boundary conditions, or combinations of these two conditions. We study in detail both the effect of the mass and the 
temperature corrections to the Casimir effect. 

We begin by computing the Casimir energy using exponential cut-off method when the massive scalar field is 
confined in a cylinder of arbitrary cross section in the background spacetime. Since a field should not have quantum 
fluctuations in the infinite mass limit [3l|, we impose the condition that the Casimir energy should vanish when 
the mass approaches infinity. We show that this condition allows us to regularize the Casimir energy and explicit 
formula of the regularized Casimir energy is given. We then consider the Casimir force in the piston setting [32| - 
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a piston dividing a closed cylinder into two chambers - which becomes fashionable nowadays. It is observed that 
the regularization procedure actually does not affect the Casimir force that acts on the piston. More precisely, the 
Casimir force acting on the piston is the same whether we compute it using the cut-off dependent Casimir energy before 
regularization or the regularized Casimir energy. The piston approach can be considered as the correct regularization 
procedure for computing Casimir force between two parallel plates 4]. From the results for piston, we deduce the 
formula for the Casimir force acting on a pair of parallel plates embedded orthogonally in an infinitely long cylinder 
[s^l- By taking the limit where the cross section of the cylinder is infinitely large, we obtain the finite temperature 
Casimir force density acting on a pair of infinite parallel plates in a (c?i + l)-dimensional macroscopic Minkowski 
spacetime, with the presence of an n-dimensional internal manifold. An advantage of our approach is that we obtain 
a formula for the Casimir force acting on a pair of parallel plates embedded in an infinitely long cylinder as a series 
over elementary functions, which enables us to derive some properties of the Casimir force easily. It is shown that the 
Casimir force is attractive if both the plates assume Dirichlct boundary conditions or both plates assume Neumann 
boundary conditions. Moreover, the magnitude of the Casimir force is always a decreasing function of the plate 
separation a and the mass m. The latter supports the assumption that quantum fluctuations vanish in the infinite 
mass limit. In the case where the surrounding cylinder assumes Neumann boundary conditions, we show that taking 
the massless limit will change the nature of the Casimir force from short range to long range. For the influence of the 
internal manifold, it is shown that the Casimir force is enhanced in the presence of extra dimensions. A stronger result 
shows that the Casimir force is an increasing function of the size of the internal manifold. By passing to the limit of 
infinite parallel plates, all these properties are preserved, although some of them is not obvious from the formulas for 
the case of infinite parallel plates. 

Besides the properties of the Casimir force, we derive explicit formulas for the asymptotic behaviors of the Casimir 
force in different limits, such as low and high temperature, small plate separation, small mass and large cross section. 
It is shown that when the plate separation is small, the Casimir force is dominated by terms that are independent 
of mass. This shows that the effect of mass is less significant if the plate separation is small. On the other hand, 
the leading order term of the Casimir force is linear in temperature in the high temperature regime. In the case that 
the size of the internal manifold is comparable to the plate separations, the behavior of the Casimir force is quite 
complicated and it depends strongly on the geometry of the internal manifold. 

Although we assume that the macroscopic spacetime is Minkowskian, the results of this paper can be easily gener- 
alized to the case where the macroscopic spacetime is also curved. By setting the size of the internal manifold to be 
zero, one can obtain the corresponding results for spacetime without extra dimensions. 

Throughout this paper, we use the units where h = c = ks ~ ^ except for the figures. 



II. BASIC FORMALISM 



We consider a background {d + l)-dimcnsional spacetime of the form 

^/di+i y. j^n^ ^Y^^^^ Afdi+i is the {di + 1)- 
dimensional Minkowski spacetime. A/"" is an n-dimensional internal space, which is assumed to be a compact connected 
manifold without boundary and d = di + n. Let the spacetime metric be given by 

ds^ ^giiudx^dx" = rjafjdx'^dx'^ - G abdy"^ dy'' , 

/X, ly —0, 1, . . . , d; a, /3 = 0, 1, . . . , di; a, 6 = 1, . . . , rt, 

where r]a/3 = diag (1, —1, . . . , — 1), j/" = for a = 1, . . . , n and Gabdy'^dy'' is a Riemannian metric on A/"". In this 

article, we are interested in studying the finite temperature Casimir effect for a massive scalar field (p{x) satisfying 
the equation of motion 

^ d^^\g>'-'d,+m^+Cm]^{x)=0, (1) 



■91 

where 91 is the scalar curvature of the background spacetime and ^ is a coupling constant — f = corresponds to 
minimal coupling and ^ = (d — l)/4(i corresponds to conformal coupling. We assume that the field (p{x) is confined in 
a cylinder of the form cyl = [0, L] x 57 x A/"", and satisfies Dirichlet or Neumann boundary conditions on the boundary 
of the cylinder. Here is a simply connected domain in M'^^^^. 

If the field ip{x) satisfies Dirichlet boundary conditions, the eigenfunctions satisfying ([T]) are given by 

ipkjA^) = e-^'^^-'' sin^^D-A^^^ ■ ■ ■ ^^"'Wy), (2) 
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where fc, j € N, ^ G N = N U {0}. The function (/>Dy (x^, . . . , x'^^) is an eigenfunction of the Laplace operator with 
Dirichlet boundary conditions on O, i.e., 



di 

4 = 2 ^ 



52 



The function '^i{y) is an eigenfunction of the operator — — ^91g = — VG ^ daVGG°-^db — i^Q^g, (9^g = is the 



scalar curvature of the metric Gabdy'^dy^) on Af" with eigenvalue ■ 



- (Ag + ^iHg) $i(2/) - ojlf.iMy)- 

The eigenfrequency w^j-^; is given by 




Since fl is simply connected, the eigenvalues ui^ j^.j are all nonzero. For Neumann boundary conditions, the eigen- 
functions are 

=e-*'^'=-'*cos^0^;,(x2,...,a;'^O'i>i(2/), (3) 

where kJJ G N = NU{0}. The function (fiN-jix"^, ■ ■ ■^x'^^) is an eigenfunction of the Laplace operator with Neumann 
boundary conditions on fl, i.e., 



d(j)N-j 



an 



where n is a unit vector perpendicular to dil. By convention, ipi\i-Q{x'^, . . . ,x'^'^) is the constant function with zero 
eigenvalue. For the eigenvalues 1^^./, we assume that w^.j > and there are exactly k of them that are equal to zero. 

The boundary conditions considered above are homogeneous boundary conditions. In the following, we are going to 
consider the Casimir force for parallel plates and also Casimir force in the piston setup. In these scenarios, we will also 
consider mixed boundary conditions. We consider two cases. For the first case, we assume that the field ip{x) satisfies 
Dirichlet boundary conditions on the wall [0, L] x dfl x A/"" and the wall x^ = 0, and Neumann boundary condition on 

the wall X — L. The eigenfunctions are obtained from ^ by replacing the function sin by sin ■ For the 

second type, we interchange the roles of Dirichlet and Neumann on the boundaries. In this case, the eigenfunctions 

are obtained from ([3]) by replacing cos by cos ^^'^^^^^ . 

Before ending this section, we define the variables R and r by i? = Yol{Qy/^'^^~^'> and r = Vo^A/")^/" which have 
the dimension of length to measure the size of the domain Q and the internal manifold A/"". Throughout this article, 
we assume that the size r of the internal manifold is smaller than any measurable length in the {di + l)-dimensional 
Minkowski spacetime. The re-scaled variables 



^Ru}n,*;j, * ^ D OT N, 

are dimensionless variables and are invariant under the re-scaling of the domain 17 and the manifold J\f" 



III. THE CASIMIR ENERGY 
A. Zero temperature Casimir energy 

At zero temperature, the Casimir energy is defined naively as the sum of zero point energies: 
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To take into account the homogeneous and mixed boundary conditions, we let 



/TY(k + a)\ 

oJkj.i =uJkj.i{a,*;m) = \ I j + ajf^^^ ,^ + uj]^.j + m'^ , k e j e J^,l e N, 

where a = 1 for homogeneous Dirichlet boundary conditions, a = for homogeneous Neumann boundary conditions 
and a = 1/2 for mixed boundary conditions; * — D (resp. * = N) when ^p{x) assume Dirichlet (resp. Neumann) 
boundary conditions on the component [0, L] x dVt x A/"" of the boundary of the cylinder; J^) = N and Jm = N. Using 
exponential cut-off regularization, we define the cut-off dependent energy as 



oo oo 



k=ojeJ, 1=0 



Using the same method as the massless case [T5l|, we find that up to the term A", the small-A expansion of Eq~^{L; A) 
is given by 

EcIs{L;\) = ^ ^^^^^^^^:jY~^Ccyi,a,*;»(m)A*"'^"i - i°^^ Ccyi,a,*;d+i(m) + ipP^^, i Ccyi.a.^s; "^), (4) 

where Ccyi.a,*:i(™) are the heat kernel coefficients defined so that 

oc oo Af— 1 

5m E e"*"'-' = E Ccyi,a,*;^(m)iT^ + O (i^) as t ^ 0+, (5) 

k=OjeJ, 1=0 i=0 

Ccyi,a,*is;m) is the zeta function 

oo oc 

k=ojeJ, 1=0 

and the finite part of a meromorphic function f{z) with at most simple pole at a point z = zq is defined by 

R-es^=2o/(2;)' 



FP.=.J{z) = lim fiz) - 

z^zo \ 

As A —> 0^, we see from (j4]) that the divergent part of the Casimir energy contains divergence of order log A and A~\ 
i = 2, . . . ,d + 1, with coefficients depending on the coefficients Ccyi,c(,*;i(TO), which can be expressed in terms of L, m 
and the geometric invariants of the manifolds f2 and A/"". In particular, the leading divergence 

r(ci + i) , _ r(rf + 1) LVoi(»)Voi(AA") 

is called the bulk divergence and is usually subtracted away in the definition of Casimir energy. The other diver- 
gences are called hypersurface divergences and regularization is required to remove these divergences. A conventional 
method, known as zeta regularization (33 |. set all the hypersurface divergences to zero and define the regularized zero 
temperature Casimir energy to be 

^cTs^'^~°(^) =^FP,^_lCcyl,a,*(s;m) -|- ^ log ^i^RcS,^. i Ccyl,a,* (s; H 
= 2^Ps=-^<^cyl,a,*(S;m) - -^-j=Ccyl,a,*;d+l[m), 

where ^ is a normalization constant with dimension length"^. This is tantamount to subtracting the divergent terms 

^-7^r7V^Ccyl,a,*;«(m)A* - ' r- Ccyl,a,*;d+l(m) (7) 

1=0 ^ y—) ^v^r 
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from the cut-ofF dependent Casimir energy ([4]) and set A = 0. However, for massive scalar fields, this definition is not 
sufficient. Since a field should not have quantum fluctuations in the limit of infinite mass [sij . it is natural to require 
the regularized Casimir energy -Eq^s vanish when the mass m approaches infinity, i.e., 



To satisfy this condition, one need to extract the leading behavior of the Casimir energy in Q or ([5]) when the mass 
is large and subtract away those terms that give nontrivial limits when the mass tends to infinity. Using the fact that 



EEE' 

k=ojeJ, i=Q 



k=ojeJ, 1=0 { 



7r(fc + a) 
I 



it is easy to show that the heat kernel coefficients Ccyi,a,*;i(TO), i — 0, f,2,... at any mass can be expressed as a 
polynomial in m~ with coefficients in terms of the heat kernel coefficients Ccyi,Q,*;ti+i(0) when m = 0. More precisely, 



Ccyl,a,*;i{m) = ^ Ccyl,a^*;i_2j (O)to^-' . 

Therefore, the divergent part of the Casimir energy ^ can be expressed as a polynomial in m^. As a result, it has to 
be subtracted away in the regularization procedure and therefore do not contribute to the regularized Casimir energy. 
For the regular part of the Casimir energy given by the second term that proportional to Ccyi.a,*-d+i{nT-) would 
not contribute to the regularized Casimir energy for the same reason above. For the term (1/2)FP^^_ i Ccyi.a,*(s; m), 
we write the zeta function Ccyi,a,*(s; to) as the Mellin transform of a heat kernel 



1 r°° °g oo 

Ccyi.a,*(s;H - — / i^-'E E E^"*""^'^^ 

V / I — n 7 i—n 



Define 



^w = EEE°^p 

k=aje.Jt 1=0 { 



■n{k + a) 



L 



k=oje,J, 1=0 



2 2 



d+l 
1=0 



Then 



K{t) =0 (t) as t^O^ 



K 



it) = o (a) 



as t — > oo. 



and we find that 



CcyLa,* (-5 j TTl) 



— t^-lg-*™' f ^ CcyLc,*;,(0)t'^ + k{t)\ dt 



r( 



d+l 

E 

i=0 



"'^cyl,Q:, + ;i(0)^ 



-2s-i-\-d 



m Jo 



k{t)dt 



gives an analytic continuation of Ccyi,Q,*(s; Tn) into the domain Re s > —1. Putting s = — f /2, we find that 

■0(-^)+logm2 _ 1 



;^FP^^_iCcyi,a,*(s;m) ^ Ccyi,a,*;<i+i(TO) +D(m) 



4\/^ Jo 



where 



^ d+l 

M = "^y^ E Ccyl,a,*;d+l-i(0)m*FP,^_ , r(s). 



Therefore, as m ^ oo. 



■0 (— i) + log 771^ 

FPs^_ 1 Ccyl,a,*(s; m) = —= Ccyl,a,*;d+l{m) + D{m) + o(l). 

2 4-^7r 



(8) 



This implies that to obtain the regularized Casimir energy that vanishes when to — > cxd, we should subtract away the 
term 

2^ ^T7i^Tr^CcyLa,*;i(TO)A' + -^-^ — Ccyl,a,*;<i+l(m) + £)(m) (9) 

i=0 ^ V 2 j ^V^r 

from the cut-ofF dependent Casimir energy ([4]) or subtract away the term 



V-l-i) +log (to//z) 
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-Ccy\,a,*-d+l{'m) + I)(to) 



from the zeta regularized Casimir energy From ([5]), we then find that the regularized Casimir energy is given by 
,,cg,T=o,,,_l^^ , ^{-\)+\ogm\ ^^^^ ^^^^^ 1 r,_3_,„. 



Ec^s {L) =-FP,^_iCcyLa,*(s;TO) - ^ c,yi,^,.,d+i{m) ~ ^{m) = J t-^e"*" K{t)dt. 

(10) 

To study the behavior of the regularized Casimir energy with respect to the variation of the plate separation L, we 
use the fact that for a = 0, 1/2, 1, 

e = - 5: e + ^ t ^e--e ^ + ^ ~ a. (11) 

k—0 k— — oo k— — oo 

This implies that 



oo 



1 \. / ^ ^ r(5-i)^ / 1 \ 



fc=l J. i=o V 
where K^{z) is the modified Bessel function of second kind, and 



CnxAf,*{s;m) = E E +'^A/-;; +™^) 

jGJ. 1=0 



As a result, we have 



Ec^f~\L) = l^FP,^_iCoxA^,*(s;™) + Rcs,=-i (r(s)CaxA^,*(s)) - FP,=_i {T{s)Cnx^^,M) 

V.(-i)+logm^ ^^^^ 1 ^^^„2....\/<-^+"^^'+™' 

*^cyl,a, + ;d+l 



(TO) - D(TO) - ^ E E E ^ , (2fcL^<,,^.+.^,+TO2 

(13) 



^ fc=i jeJ. i=o 



Using pT|) . we also have 



00 00 



E E E e-*"^- - (^^- + J - " + -d. ) E E e-K.^^--'-"^) 



(^^t-^ + i - a + c.d.^ 1^ E cny.H.*-A'm)t~^ + O {t^^-^^j as i ^ 0^ 
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where e.d. is the exponentially decay terms, and ca-^N' are heat kernel coefScients for an elliptic operator on 
r2 X A/". This implies that 

Therefore, we read from (|13p that the regularized Casimir energy can be written as a linear polynomial in L, plus a 
term that decays exponentially as L approaches infinity, i.e., 



(L) = So + Sli - ^ E E E e^"''" ^ ^1 (2fcLy'c.,V,^.+c.^,+n.2) . (15) 



B. Finite temperature Casimir energy 

Now we take into account the temperature correction to the Casimir energy given by 

T^log^l-e-"'--^-'/^) . (16) 

This summation is finite and no regularization is required. Therefore, the cut-off dependent finite temperature Casimir 
energy is defined as 



It is well known that this can be computed in terms of the finite temperature zeta function given by 



oc oo 



fe=ojeJ, 1=0 p=-oo 



More precisely (see |15l|). 



p (r u f , log[A/i]-V'(l)-log2 + l 

Ecasm^) T^fd-i\ Ccyl,Q,*;j(m)A H Ccyl,Q.*;rf+l(TO) 

1=0 ^ y—) ^^'^ 

T 

- (Ccyl,a,*;T(0;TO) +log[Ai^]Ccyl,a,*;T(0;m)) . 



(17) 



For the regularization, we observe that the temperature correction to the Casimir energy (jl6p vanishes when the mass 
TO approach infinity. Therefore, the large-m leading behavior of the cut-off dependent finite temperature Casimir 
energy (jl7p is the same as the large-m leading behavior of the cut-off dependent zero temperature Casimir energy 
given by ([9]). After subtracting away these large-m non- vanishing terms, we find that the regularized finite temperature 
Casimir energy is given by 

y log[^jVv(l) 

£^cTs(i) = - 2 (CcyLa,*;T(0;™) + log[/i']Ccyl,a,*;T(0; m)) Ccyl,„,*;rf+1 (m) - S)(to). ^^^^ 

As a remark, in the zeta regularization scheme, the regularized finite temperature Casimir energy is defined as 

^Zlm = -\ (Cyl,a,*;T(0;™) + l0g[M']Ccyl,o,*;T(0;m)) . (19) 

This expression does not go to zero as the mass to approaches infinity. However, if we subtract away from (fTO|) the 
large-m leading terms that give nontrivial limits, we would obtain p^ . We would like to remark that although the 
expression (jlSp for the regularized Casimir energy is derived for the case where the field is confined within a cylindrical 
cavity, this formula is in fact valid for an arbitrary closed cavity. Similarly, the formula for the regularized Casimir 
energy at zero temperature given by (fTU]) is also valid for any closed cavity. We would also like to point out that 
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the high temperature expansion of the zeta regularized finite ternperature Casimir energy for massive scalar field in 
manifolds with boundaries have been obtained in [s^ [s^, [3 7| . [s sl. [soj . 

Using the same method as in the derivation of (|T2|) (see [l5l|)" we find that 

^ s Ccyl,a,*-d+lim) 
Ccyl,a,*;T(0;m) = , 



oo oo 



27rika 



Ccyi,a,*;T(0; ™) =Ao + AiL + ^ ^ ^ J] -j— exp (^-2kL ^ iuf, ,^ ^ + + m2 + (2^pr)2) 



k=l 1=0 p=-oo 

where Aq and Ai are independent of L, and Ai is given by 



Ai = - ■7^cnx7V,*:d+i("i) + 7^FPs^_i {r(s)Cj2xA/-(s; ™)} 



2 oo oc 



-^-l^l^l^ ^1 T 

]£J, 1=0 p=l ^ \ 

Together with (fT^ , we find that the regularized Casimir energy can be expressed as the sum of a linear polynomial 
in L plus a term that decays exponentially as L approaches infinity, i.e., 



Ec!siL) = So + Si^ - 2 E E E E (-2fci V-'^l*., + ^'m, +^' + i^^pTy) , (20) 



where 



1 ''^^ 1 
Si ^-^V" FP^^_ir(s) cnxM,*:d+i-i{0)m^ - t— {FPs=-i {r(s)CoxA/'(s; to)} + log[TO2]cf2x7V,*:rf+i(TO)} 
an ^ — ' L 2 J j^Ti" 

4=0 

T ^ ^ ^ ^J^k*-j+^lr;l+^^ ( P\l^fx*-,3+^M;l+'^^^ 

jSJ. !=o p=l ^ \ 

Before ending this section, we would like to comment that it is easy to deduce from the results above that the 
Casimir energy for mixed boundary conditions (a = 1/2)) are related to the Casimir energy for homogeneous boundary 
conditions (a = or 1) by 

i^cas [L-a^\,^ = D\ =Scas (2i; a = 1, * = Z?) - E^^, {L;a^l,* = D), 

Ecs [L;a= N] -Scas (2i; a^O,* = N)- Ec^, {L;a^O,* = N). 



C. Special Case 

Here we consider the special case where the cross section of the cylinder is a rectangular region [0, L2] x . . . x [0, L^jJ 
and the internal manifold Af" is an n-torus T" - a product of n circles with radius ri , . . . , r„ respectively. Notice that 
the scalar curvature of the torus T" is zero. The spectrum uj^ ^.j of the Laplace operator with Dirichlet or Neumann 
boundary conditions on f2 is given by 

and the spectrum uj'j^.i of the Laplace operator on T" is given by 



ri I \rn 
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For simplicity, we only consider the cases where a = or 1, i.e., all the walls of the rectangular cavity [0, Li] x . . . x 
[Ojidi] assume Dirichlet boundary conditions or all the walls assume Neumann boundary conditions. The results for 
the case where a = 1/2 can be obtained from the results for a = or 1 by ((22|) . For the finite temperature Casimir 
energy inside the cavity [0, Li] x [0, L2] x . . . x [0, L^J x T" given by p?)) . we have 



cy\, D/N;i[^) — '71 Ccyl^D/N;i-2jt 



(23) 



where 



Ccy\, D/N;i 



7 - 



.1=1 J 



Sdi- 



0<i<di, 



and 



E 



Lai ■ ■ ■ La 



l<ai<...<aj<di 



is the (hyper)-surface area of j-dimensional hyperplanes. For z > di + 1, Ccyi^D/N;i — 0- The finite temperature zeta 
function Ccyi,_D/Ar;T('7i) Can be written as a sum of inhomogeneous Epstein zeta function, i.e.. 



l<ai<...<ai<di 



TT 1 1 

, — ,27rr;m , 



1 ■ ■ ■ 1 7- 1 1 



Lai 1"1 



where 



Zj{s;ai,...,aj;m) ^ ^ 



The properties of the inhomogeneous Epstein zeta function Zj (s; ci, . . . , Cj; m) has been discussed quite extensively 
(see e.g. [IS El IS Ell)- In particular, we obtain from [il] that if j is even, 

(—1)5 nim^ 



Zj{0; fli, . . . , aj; to) 
Zj(0;ai, . . .,aj]m) 



(i)! 



(-1)5 nim^ / , , , /?■ + 2 , 

' logTO^ + V' ( ) -ip{l) 



(§)! 



2m 2 



E E 

(fei,...,fcj)6Z3\{0} \i=l 



2\ 4 



Ki 27rTO, 



(24) 



and if j is odd, 

Zj(0; ai, . . . , a^; to) =0, 



Z (0; ai, . . . , Oj; to) =t — ■ tF 



2to2 



E E 



2\ -1 



27rTO, 



(25) 



Notice that the heat kernel coefficients (j23p are non- vanishing when m — > 00. Therefore, they are subtracted away to 
obtain a regularized Casimir energy that vanishes when m — > 00. On the other hand, we obtain from (j24p and 
that 

-| (Ccyl,D/iV;T(0; + log[M']Ccyl,D/iV;T(0; to)) - *P + i?Ca.;D/JV (^1 ) ' 
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Left 
Chamber 





Right 
Chamber 



L-a 



FIG. 1: A movable piston inside a closed cylinder divides the cylinder into two chambers. 



where ^ is the large-m non- vanishing term and E^q^^.£i/j^{Li) is the regularized Casimir energy given by 



E. 



Cas-D/N^ 



E 

l<(Ti<...<o-i<di 

+ n + l 



E 



(A:^,,...,fc^^,/i,...,/^,p)eZ^+- + i\{0} 



n 

X I J2 [i-.fc-J' + Et'^'^j^j 
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(26) 



For massive scalar field with Neumann boundary conditions, the regularized Casimir energy is negative. 

The terms with p = in (|26p give the zero temperature Casimir energy. The sum of the terms with cti > 2 
corresponds to the term Eq in (PO)) which is independent of Li. The sum of the terms with <Ji = 1 and fco-i = fci = 
corresponds to the term Si in ((20)l which is proportional to Li. The rest of the terms decay exponentially when 
Li — > cxD. By applying the formula 



£ exp (^-^ - t^^^ dt - 2 (^^y X,(2/37) 



exp {—tk"^) = \fT^t 2 ^ exp ( — 



and the formula 



t 

k— — oo k— — oc 

we can rewrite the regularized Casimir energy £^^^^g(L) ([26| in the compact form ([201 



(27) 



IV. THE CASIMIR FORCE 

A. The Casimir force acting on a piston embedded in a closed cylinder or on two parallel plates embedded 

in an infinitely long cylinder 



To find the Casimir force acting on the walls of the cylinder, one need to take into account the Casimir energy of 
the outside region, which is not easy to compute. However, there is a setup closely related to the cylinder where the 
Casimir force can be computed. As in [33 |. one consider a freely moving piston (see FIG.[T]) dividing the cylinder into 
two chambers: the left chamber [0,a] x 17 x A^" and the right chamber [a,L] x il x A/"". In this case, the Casimir 
energy of the region outside the cylinder does not have any effect on the piston. The Casimir force acting on the 
piston is the sum of the Casimir forces arise from the variations of the Casimir energies in the left and right chambers: 



Fcas{a;L) 



d_ 

da 



(EcUa) + EcUL - a)) ^ 



(28) 
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From (|20p . we find that the contribution to the (regularized) Casimir force from the left chamber is given by 



jeJ. /=o p=-oo exp {2a^ul^^^^ + + + {2npT)^j - (-1)2" 

where Ei is a term independent of a given by (HJ). The contribution to the (regularized) Casimir force from the 
right chamber is negative the contribution from the left chamber, with a replaced by L — a. As a result, the term Ei 
cancels out and we find that the Casimir force acting on the piston is given by 

FcUa; L) = ^^Ss(a) " PgLsiL ~ a), (30) 

where 



^c°:.(«) = -(-if"7^EE E — 7 — / N (31) 

jeJ. 1=0 P=-oo exp [2a^ujl , .j + uj]^j + + {2npTfj - (-1)2" 

can be interpreted as the limit of the Casimir force when the right chamber is infinitely long. 

Since the difference between the cut-off dependent Casimir energy (jl7p and the regularized Casimir energy (j20[) is 
given by ©, which by is a linear function in L, this implies that we will obtain the same result for the Casimir 
force whether we use the Casimir energy before or after regularization. In other words, the Casimir force acting on 
the piston is independent of the regularization procedure employed. 

It is easy to deduce from (1311) that the Casimir force for a = 1/2 is related to the Casimir force for a = or 1 by 

^^Cas [a- L; a = i; * = D/N^ = 2Fc,s (2a; 2i; a = 1/0, * = D/N) - F^s (a; L;a = 1/0, * = D/N) , (32) 
which also follows from ((22|) . 

The expression ([5T|) is negative if a = or 1 and positive if a = 1/2. Moreover, its absolute value is a monotonically 
decreasing function of a. Therefore, if the quantum field assumes either Dirichlet or Neumann boundary conditions on 
both the piston and the walls of the cylinder, the Casimir force acting on the piston is an attractive force tending to 
pull the piston towards the closer wall. If the quantum field assumes Dirichlet (resp. Neumann) boundary conditions 
on the piston and Neumann (resp. Dirichlet) boundary conditions on the walls of the cylinder, then the Casimir force 
acting on the piston is a repulsive force tending to restore the piston to its equilibrium position = L/2. In both 
cases, the magnitude of the Casimir force increases as the piston is moving away from the equilibrium position. 

Eq. (j3ip can also be interpreted as the Casimir force acting on two parallel plates embedded in an infinitely long 
cylinder with cross section x A/"". It shows that the Casimir force between the plates is attractive if the field 
assumes the same boundary conditions on the plates and is repulsive if the the field assumes different boundary 
conditions on the plates, regardless of the boundary conditions assumed on the surrounding transversal wall. Writing 
<^n,*y = ^.j/R and uJ^f■l = w^.j/r, eq. (pij) shows that when a/r or a/R or am is large, the magnitude of the 
Casimir force decays exponentially. At high temperature T, eq. ([3T|) shows that the leading term of the Casimir force 
is given by a term linear in T: 



exp [2a^ujl^^ + ujj^ i + m^j - (-1)2" 

and the remaining term decays exponentially. Notice that if we consider the contribution to the Casimir force from 
the left chamber , then the result (jAip of Appendix |^ shows that when T 1 , 

« '"'*(«) (^) Cr{^ + 2)caxA^,M-i-.(™)(2T)^+2 + 0(T). (34) 

We observe that in general, there are terms of order T^, T^, . . . , r^+^. In particular, the leading term is 

TT 71-2 
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which is the Stefan-Boltzmann term. Usually this term is subtracted away in the regularization of the Casimir energy 
since it can be interpreted as the contribution to the vacuum energy in the absence of boundaries. However, besides 
the limiting case of infinite parallel plates, there are still terms with order T^, . . . , T"^. Nevertheless, since these terms 
are independent of a, they cancel the corresponding contributions from the right chamber and the high temperature 
leading term of the Casimir force acting on the piston is a term linear in T given by (j33[) . This is usually called the 
classical term [i^ |4^, |43, 13 due to the absence of the Planck constant h in this term. In a recent work on Casimir 
effect of electromagnetic field in three dimensional ideal metal rectangular box fW] , it has been argued that the terms 
of order T^, . . . , T'^^^ have to be subtracted away in order to be consistent with thermodynamics. Here we find that 
in the piston scenario, such terms are naturally absent due to the cancelation between the two sides of the plate. 
For the low temperature asymptotic expansion of the Casimir force, we use the formula 



exp 
1 d 



(2ay'<,^^.+a;2,_, + m2 + (27rpT)2] - (-l)^^ 



(35) 



2-^7? da I 
and the formula (|27|) . These give 
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In the zero temperature limit, we obtain the zero temperature Casimir force from the terms with p = 

^ oo oo 

V)4EEE 
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(37) 



which can also be derived directly from eq. (jlSp . The temperature correction to the Casimir force is the sum of the 
terms with p ?^ in p6)) . It goes to zero exponentially fast when T —> 0. Using the identity -fC2(z) ^ Ko{z) + 2Ki{z)/ z, 
we can rewrite the zero temperature Casimir force ([37]) as 



^CaL («) 



^EEE^^^"" 

k=ijeJ, 1=0 




Ki [2kaJu;f^ 



{uJn,*;j + ^lr,i + m"^) Ko (2A:ayC^~ 




(38) 



which shows manifestly that the zero temperature Casimir force is attractive if a = or 1 and is repulsive if a = 1/2. 
For the behavior of the Casimir force with respect to the variation of mass, we observe that the function 



is a decreasing function. Therefore the magnitude of the Casimir force acting between a pair of parallel plates 
embedded in an infinitely long cylinder (of arbitrarcy cross section) is a decreasing function of m when a = or 1. In 
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other words, for a pair of parallel plates with identical boundary conditions, the increase in mass reduces the Casimir 
effect. This property cannot be obviously inferred from the expressions for the zero temperature Casimir force given 
by (|37p or ([38|l . We see here the advantage of considering the Casimir force at any finite temperature. It enables us 
to derive some properties of the Casimir force from the expression (1311) . which cannot be derived directly from the 
expression for zero temperature Casimir force (jST]) or (|38|) . 

For a pair of parallel plates with mixed boundary conditions or when a — 1/2, the situation is different. The 
function 



1 



is increasing when x e [0, 1.2785] and decreasing for x > 1.2785. Therefore, when the mass increases, the magnitude 
of the Casimir force may first increase and then decrease exponentially. 

Notice that when the mass m decreases to zero, ([5T|) naively shows that the Casimir force tends to the Casimir 
force for massless scalar fields [lH| • To be more careful, we need to discuss the cases where the field assumes Dirichlet 
boundary conditions and Neumann boundary conditions on the surrounding walls [0, L] x x A/"" separately. If the 
field assumes Dirichlet boundary conditions on the wall [0, L] x dfl x A/"", then the Dirichlet eigenvalues ujq jj.j are 
all nonzero. In this case, we can immediately set m = in PT|) and obtain the Casimir force for massless scalar field. 
However, if the field assumes Neumann boundary conditions on the wall [0, L] x dft x A/"", there is exactly one zero 



Neumann eigenvalue ^.q. 



In this case, there are k pairs of (j, I) such that cj^ j^.j + oj'j^.i — 0. Separating the sum 



over (j, I) e JAf X N in ((31 
we find that the we can immediately set 



as a sum over those (j, I) with uf- 



m 



^o^wu + ^Ifa and those {j, I) with w^^jvy- + ^^If-i ^ 0, 
in the latter sum to obtain the contribution of the modes with 
^n.N-j + 7^ to the Casimir force due to massless scalar field. For the sum over (j, I) with i^n^N-j + ^"If-.i 
given by 
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P— — 00 



exp (2ay/m^ + (27rpT)2^ - (-1)^" 



(39) 



we have to be careful when taking the massless limit for the term with p = 0. Using the fact that 

X „, . X 



l + 0(x). 
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we find that the massless limit of ([55)1 is given by 
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as X ^ Q, 
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(40) 



Except for the factor k, this is the finite temperature Casimir force between a pair of parallel plates in (1+1)- 
dimensional Minkowski spacetime due to massless scalar field with Dirichlet or Neumann boundary conditions on 
both plates (for a = or 1) and with Dirichlet boundary condition on one plate and Neumann boundary condition 



on the other plate (for a = 1/2) (see [50|). 
(|40ll has an alternative expression given by 



As in 



the Casimir force due to the modes with jy.^- 
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(41) 



and 
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(42) 
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The first term in (|¥T|) and (|1^ give the respective zero temperature Casimir force. Notice that in contrast to the 
massive case where the force decays exponentially, they decay in the order 1/a^ when a is large. This is a long 
range force which is the subject of study in the context of electromagnetic fields in the recent work [5l|. At finite 
temperature, (|40p shows that this long range force is present if and only if a ^ 1/2, i.e., if and only if the boundary 
conditions assumed on the two plates and the walls of the cylinder are both Neumann conditions. In this case, the 
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long range force is of order T/a. From this analysis, we find that long range force may exist only in the massless 
case when Neumann boundary conditions are assumed on the wall of the cylinder and the two plates or when the 
temperature is zero and Dirichlet boundary conditions are assumed on the two plates. The transition from massless 
to massive field will change the nature of the force from long range to short range. 

To investigate the dependence of the Casimir force FQ^^{a) on the size R of the cross section f2, we rewrite 
as Lo'^ where the re-scaled frequency ^.j is independent of the size R of the cross section il. Using the same 

argument about the dependence of the Casimir force on mass, we see that if a = or 1, the Casimir force increases 
when the size R increases. The asymptotic behaviors of the Casimir force when the plate separation a is much smaller 
than the size R of the cross section 57 is derived in Appendix [B] We read from (IB2[) and (IBSP that if a ^ i?, then if 
aT> 1, 
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(43) 



and if aT < 1, 
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(44) 

Here CQju ^f-i are the heat kernel coefficients of the Laplace operator with Dirichlet (* = D) or Neumann (* — N) 
boundary conditions on Q./R. Notice that when r a R, the large-i? non- vanishing terms of the Casimir force 
can be written as a polynomial of order di — 1 in i? with coefficients depending on the geometric invariants co^,;^ of 
il, and Bessel series that depend on the geometry of the internal manifold A/"", and the plate separation a. From 
these expressions, it is easy to read that if the mass m is also very large, the Casimir force decays exponentially. In 
the case the mass m is small, or more precisely if am <C 1 ^ Rm, we obtain from (jB4j and (|B6|) that if a = or 1, 
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/TTQ '■ — ' V a / '■ — ' q. 
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(46) 



if T = 0. The omitted terms goes to zero when am ^ or R/a — > oo. Notice that these asymptotic behaviors contain 
logarithmic terms in am which goes to zero when the mass m approaches zero. In the massless case, we have 

F^^^ia; m^O)^ ^{d, - z - l)r Ud, - i) - ^'Jhl^ + o(«0)^ as r « a - 0+ 



if aT > 1; and 

Fc.sia;m = 0)^~^'f:^id,-^)T(^^^^^ as r « a ^ 0+ 

if aT ^ 1. From ([45]) and (|46|) . we find that when m is small, then as a ^ 0+, the leading order term of the Casimir 
force is given by 

F^^M - ~{d, - l)r (^) Cn{d,)^^^, if aT » 1, 

V ^ / (47r) 2 a 

and 

FSL.{a) ^ -d,r Ud, + 1);^^^, if aT « 1, 

\ z J (47r)T~ 

respectively, which are independent of the mass m. This shows that when am <^ 1, the mass correction to the Casimir 
force is not significant. 

If rfi = 3 and a = or 1, ([iS]) and give respectively 

as r ^ a ^ 0+ 

(47) 
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(48) 
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if T <C 1. If a = 1/2, similar computation gives 



m2 log(am)2 + ''^/f*'^^ \og{amf \ + 0(a°) 



as r <C a 




Ci?(3)H 24^2 4^y7 log(am) S + 0(a") as r < a ^ 0^ 
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OTT lOTT 

if r ^ 1. When di = 3 and fl = 0, Lg] x [0, L3] is a rectangle, (|48p gives the correct behavior of the Casimir force 
when a <C 1 which was derived in 4J] . 

In (j46[) . we only give the asymptotic behavior of the zero temperature Casimir force when am <C 1 ^ Rm. At 
finite low temperature, the behavior of the thermal correction depends on the relative strength of the temperature T 
and mass m. If T ^ to, the thermal correction term can be expanded in a power series of to^ plus an exponentially 
suppressed term. If T 3> to, then the thermal correction term can be written as a sum of a power series in as in 
the case of T <C to and the term (|B7[) . which contains logarithmic terms of m/T. 

For the influence of the extra dimensions, we rewrite lj^-i as w^.j/r, where the re-scaled frequency i^^.; is inde- 
pendent of the size r of the extra dimensions A/"". We find from each of the formulas and asymptotic expansions 
for the Casimir force F^^{a) derived above that, as the size r is very small compared to the plate separation a, the 
terms with nonzero tf^.j contribute Casimir force that are exponentially small. In the limit the internal manifold 
A/"" vanishes, i.e. r ^ 0^, only the k terms corresponding to ll>]^.i = remain and they give k times the Casimir 
force in (di -|- l)-dimcnsional Minkowski spacetime. In other words, if k = 0, the Casimir force goes to zero in the 
limit of vanishing internal space. This is definitely not a desired physical situation. The situation that is of physical 
interest is the recovery of the Casimir force in (di + l)-dimensional Minskowski spacetime in the limit of vanishing 
extra dimensions, or equivalently k = 1. This happens in particular when A/"" has zero scalar curvature or when ^ = 
(minimal coupling) and A/""" is connected. In this case, (ISTI) shows that the presence of extra dimensions enhances the 
Casimir effect. In case a = or 1, the expression PT|) shows that the magnitude of the Casimir force becomes larger 
when the size of the internal manifold r is larger. 

When the size r of the internal manifold A/"" becomes comparable to the plate separation a, the correction to the 
Casimir force in (di + l)-dimensional Minskowki spacetime due to the extra dimensions (i.e. the contributions from 
the terms with uj]^.i ^ 0) can become substantial and it depends on the geometry of the internal manifold A/"". 



B. The Casimir force density acting on a pair of parallel plates 

In this section, we consider the limit where the size R of the cross section fl goes to infinity, which is tantamount 
to two parallel plates in a (di -I- l)-dimensional Minkowski spacetime M'^^'^^ , with an n-dimensional internal manifold 

compactified to Af" at every point of M'^^^^. In this case, we should consider the Casimir force density J^Q^^^ia) on 
the plates — and = a which is defined as the limit 

From (|43|) and ((44|) . we observe that when R is large, the leading term of -Fcas(a) is of order R'^^~^ coming from the 
term with i = 0. All the remaining terms are of order smaller than R'''^~^ . Using the fact that cn/_R,*;o — l/(2\/7i')'*^~^i 
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we obtain immediately the following high and low temperature expansions for the Casimir force density ^cas(^)' 
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(49) 

The first two terms in (|49p give the zero temperature Casimir force density and it agrees with the result obtained 

in [16]. Notice that since the Casimir force density ^casi'^) derived as a limit of the Casimir force J^^g(a), it 
follows that for a pair of infinite parallel plates with Dirichlet or Neumann boundary conditions on both plates, the 
Casimir force is attractive. For a pair of infinite parallel plates with Dirichlet boundary condition on one plate and 
Neumann boundary condition on the other plate, the Casimir force is repulsive. For either homogeneous or mixed 
boundary conditions, the magnitude of the Casimir force is a decreasing function of the plate separation a, but it 
is enhanced by the presence of the extra dimensions. If both the plates assume the same (Dirichlet or Neumann) 
boundary conditions, the Casimir force density ^Qi^s{a) is also a decreasing function of mass. In the high temperature 
regime, the leading term of the Casimir force is linear in temperature. In the low temperature regime, the Casimir 
effect is dominated by the zero temperature term. As is shown in Appendix iBl the behavior of the thermal correction 
is rather complicated and it depends on the relative strength of T and m. If T to, the thermal correction term to 
the Casimir force density is exponentially suppressed. However, if to <C T <C 1, we find from (jB7p that the leading 
behavior of the thermal correction term is 



V^Tto^^ f-log (^)\^( ^) - V.(l)) - f-^ ) TO^-^ \, if dr is even; 



(i|L)! ^ V °\T/ ^ V 2 J ' ' 7 2 



. (-^) TTO^^ . (log {£,y ~ , (^) - .(!)) |, . is Odd. 

In the limit when the mass to goes to zero, this shows that the leading order term of the temperature correction to 
the Casimir force density acting on a pair of parallel plates is 

r(^)cR{di + i)T'^^+\ (50) 



di + l 
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This is in contrast to the case where the area of the cross section is finite, in which case the thermal correction term 
to the Casimir force is exponentially small when the temperature is low. 



C. Special Case 



Here we consider the special case as in section UlI C[ where the cross section of the cylinder 57 is a rectangular region 
[0, L2] X ... X [0, Ldi] and the internal manifold A/"" is an n-torus T" - a product of n circles with radius ri, . . . , r„ 
respectively. For simplicity, we only consider the Casimir force F^^^{a) acting on a rectangular piston moving freely 
inside a semi-infinite long rectangular box when both the piston and the walls of the rectangular box assume the same 
(Dirichlet or Neumann) boundary conditions. The result for the case where they assume different boundary conditions 
can be obtained using (f32|) . From ([26]) and (I29p . we find that the left closed chamber contributes the Casimir force 
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(51) 



For Neumann boundary conditions, it is obvious that Si is negative. Using (|27p . we can rewrite Ei in the form (|2ip . 
where the terms with p = in (jSip are the T- independent terms correspond to the first two terms in \21\ . As is 
shown in AppendixEl when T is large, Ei is negative and dominated by a term proportional T'^^^ . Therefore, when 
T is large, there is a large repulsive force due to the vacuum fluctuations of the scalar fleld in the left chamber that 
tend to push the piston away from the closed end of the rectangular box. However, the vacuum fluctuations of the 
field in the right chamber give rise to a larger force in the opposite direction. The sum of these two forces is 
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(52) 
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which tends to move the piston towards the closed end of the rectangular box. The expression ([S^ shows that the 
Casimir force F^l^la; D/N) is exponentially small when any of the parameters a, m, T is large or any of the parameters 
L2, ■ ■ ■ , , ri, . . . , r„ is small. There are a few alternative expressions for the Casimir force F^^{a; D/N) which can 
be used to study the behaviors of the Casimir force at other limits. Using the formula ([M)) . we find that if a and m 
are large and L2 , . . . , L^i . . . ,rn,T are small, the Casimir force can be computed using the formula 



T7100 



-.00 oc 

(«;^/^) = iE E EE 



fc=i 



1=2 



Li 



2k^a^ 




19 



In the zero temperature limit, we find that the zero temperature Casimir force is given by 
Fgi,{a;T^O;D/N) 
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This expression shows that the zero temperature Casimir force is exponentiahy small when any of the parameters 
a, m is large or any of the parameters L2, . . . , Ld^ , ri, . . . , r„ is small. For the situation that we are more interested 
in, i.e., the case where <C a <C Lj, 1 < « < n, 2 < j < di, the Casimir force F^^{a;D/N) can be written as the 
sum of two terms, where the first term 
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is the dominating term, and the second term 
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decays exponentially when the area of the cross section L2 ■ ■ ■ Ld^ is large. The term Si in ((53l) is equal to the sum 



s^= E 

2<(Ti<...<o-i<di 



The expressions (j53p and (|54p are suitable for investigating the high temperature behavior. The high temperature 
leading term is linear in T, which is equal to the sum of the terms with p = 0. At low temperature, we have the 
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following alternative expressions for ((53)) and ([54]) : 
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The p = in (j55p and (j56p give the zero temperature Casimir force. In the limit the cross section is large, i.e. Li oo 
for i = 2, . . . , c?i, we find h'om the term with i = di — 1 in (j53p and (j55p that the Casimir force density acting on a 
pair of infinite parallel plates is given by 
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FIG. 2: The Casimir force Fc^is(a) as a function oi r/a for different values of m and T. 
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FIG. 3: The Casimir force -Fc£s(a) as a function of temperature for different values oi r/a and ' 




In Figures [2l [3] and |4l we show graphically the behavior of the Casimir force (a) acting on two parallel plates 
embedded in an infinitely long rectangular cylinder in a (4 + l)-dimensional background spacetime with an extra 
dimension curled up to a circle of radius r. We assume that the massive scalar field is subjected to Dirichlet boundary 
conditions on the walls of the cylinder and the two plates. Figure [5] shows the variations of the Casimir force as a 
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mass (eV/c^) 

FIG. 4: The Casimir force Fc^^s(a) as a function of mass for different values of r/a and T. 



function of the size r of the extra dimension, or more precisely, the ratio of the size r to the plate separation a. It 
shows that the magnitude of the Casimir force is increased if the size of the extra dimension is increased. When 
r/a > 0.4, the existence of extra dimension can contribute substantially to the Casimir force. Figure [3] shows that 
the magnitude of the Casimir force is an increasing function of temperature. For a = 6/^m, the graph shows that the 
Casimir force depends linearly on temperature when T > 2QQK. Figure [4] shows that the magnitude of the Casimir 
force decreases when the mass m increases. 



V. CONCLUSION 



In this paper, we have investigated the Casimir effect for massive scalar field with general curvature coupling in 
{d+ l)-dimensional spacetime with n — d~di extra dimensions. We consider the cases that the field assumes Dirichlet 
or Neumann boundary conditions on two parallel plates embedded in an infinitely long cylinder. We derive a general 
expression for the regularized Casimir energy that vanish in the infinite mass limit. A lots of the properties of the 
Casimir force are similar to the massless case. In particular, if the field assumes Dirichlet or Neumann boundary 
conditions on both plates, then the Casimir force acting on the plates are attractive. If the field assumes Dirichlet 
boundary condition on one plate and Neumann boundary condition on the other plate, then the Casimir force is 
repulsive. Passing from massless to massive, we find that the strength of the Casimir force is reduced if both plates 
assume the same boundary conditions. 

For the influence of the extra dimensions, we find that the presence of extra dimensions enhances the Casimir effect. 
When the size of the internal manifold shrinks to zero, one obtains the Casimir force in the {di + l)-dimensional 
spacetime if an only if a certain elliptic operator on the internal manifold has exactly one zero mode. This happens 
in particular in the case of minimal coupling or zero scalar curvature. 

For the extension of this work, it would be interesting to consider the general case of Robin boundary conditions as 
in [T6| . This will be important for considering the finite temperature Casimir effect in the braneworld model especially 
for the radion field stabilization mechanism, as has been shown in [l6| for the zero temperature case. The work along 
this direction will be reported elsewhere. 
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APPENDIX A: HIGH TEMPERATURE ASYMPTOTIC BEHAVIOR OF Ei 



In this section, we consider the asymptotic behavior of Si ((2T|) at high temperature. The first two terms in ((2T|) 
are independent of T . For the third term, we have 
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APPENDIX B: THE ASYMPTOTIC BEHAVIOR OF THE CASIMIR FORCE Fg'^^{a) IN DIFFERENT 

LIMITS 

In this section, we derive the asymptotic behavior of the Casimir force F^g(a) in different hmits. We only consider 
the case where a = or 1. The case a = 1/2 can be derived analogously. 

First we consider the asymptotic behavior when r < a < R. In the high temperature regime, it follows from (j35p 
that 
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where is the zeta function 
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and co^,;,; is the heat kernel coefficient of the Laplace operator with Dirichlet/Ncumann boundary conditions on Q, 
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This expansion only gives the behavior when am 2> 1. For the expansion when am <C 1, we go back to (jBip and 
consider the expansion of 

Tii;a) = £^aj^ t^-'J2J2 E ^ tfcVj d^j (B3) 

when am ^ 1. Recall that there are k zero eigenvalues t^^^./- We write (|B3p as the sum of Ti(i; a) and T2(«; a), where 
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For the term T2(i; a), we can use the Taylor expansion of e ™ to write 
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For a), notice that when s is large enough, and am <ti 1, 
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By analytic continuation, this formula holds for all s. If s > 1, we can put in directly the value of s in each of the 
terms. For s < 1, we have for j = 0, 1, 2, . . ., 
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Therefore, when r < a < R and am ^ 1 ^ i?m, the Casimir force F^^{a) has the asymptotic expansion 
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In the low temperature regime, we use (jlip to transform p6p to 



00 00 00 f 

EEEexpf 

2v^i ^-EEEE' 

•^0 fc=i jeJ, 1=0 p=i 

^ poc 00 00 ' ■ 



47r 



'0 



E EE^^^p' 





Trfe 




a 




+ TO 





- tk'^a^jdt 
V2T/ 



Using the same method as we derive (jB2[) . we find that 

00 

^^Ss(a)'-E^^/fl.*;*^''"'" 



'i=0 



47r 9a 



/.oo 00 00 



exp ■ 



• dt. 



- tk^a^ \ dt 



fe=l (=0 



00 00 C30 



■2V^/ i^EEE 



exp 



Trfc 



a 



47r 



f di — 1 — 



26 



i=0 ^ ' 



di — 1 — i f -I oo oo 



27r 



EE 



fe=l i=0 



((ii - i) 



d-i + l- 



+ 2 



oc oc oo 



2^.i(aT)^^^^ft^ 
;s=i i=o p=l 



-^(7rfc)2 + {au>j^.iY + [amy 



di -«-2 



(B5) 



K 



di-i-2 



+ 



(2aT) 



^1 + OO OO 



27r 



■EE 



^r— \Ta 

di-i + l 



1=0 p=l 

In case T = and am -C 1 <C -Km, 



FrCX>,J=U/' \ 



[di + il 



j=0 



2j-2 



1 (-1)^ 



a! 



X (2g - 2j + l)n^J-''"-'T (^-j + l+qj Cr{-2j + 1 + 2«) | + ^ ca/R,..d^-2j 
~ ~ I) + ^E ^(«m)2'^(2g - 2i)7r2i-2«r (-J + g) U-2j + 2q)\ 



2i-i 



47ra2 



E 

i=di+2 



g=0 ^" 



2 (am)' 



TT 



Cr {i-di + 2q) 



Tcfi E 



27ro2 



i? 

Ca,*;i I — 



i=0 



di — 1— 2 oc 



9=0 



dl + 3-i-2g 



A; 2 2 I 



(B6) 



For the temperature correction terms, it can also be written as a power series in if T < m. However, if m <C T, 
we have a contribution 
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coming from the terms with (^"j^.i = in the last term of (|B5|) . It contains logarithmic terms of m/T. 
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